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THE  FORMULATION  OF  SELECTED  FREE  BOUNDARY  PROBLEMS 
AS  CONSERVATION  LAWS 

Joel  C.  W.  Rogers 
Laurel,  Maryland,  U.S.A. 

1.  Introduction  and  Examples 

In  this  brief  paper,  we  shall  summarize  some  recent  work  on  free 
boundary  problems  which  may  also  be  considered  as  conservation  laws. 

This  Is  by  no  means  a  review  article,  and  at  best  It  Is  a  survey  of  cur¬ 
rents  of  thought  which  have  emerged  In  work  In  which  the  author  has  par¬ 
ticipated,  mostly  In  collaboration  with  others.  This  may  explain  the 
rather  parochial  nature  of  the  list  of  references  at  the  end.  The  present 
state  of  our  knowledge  Is  quite  rudimentary,  and  It  will  be  apparent  to 
the  reader  that  conjectures  considerably  outnumber  theorems.  The  subject 
matter  of  this  paper  will  be  fairly  broad;  the  price  of  this  breadth  will 
be  a  lack  of  depth. 

Generally  speaking,  a  set  of  differential  equations  will  be  said  to 
be  In  conservation  form  If,  by  formal  Integration  of  the  equations  over  a 
region,  one  can  obtain  equations  for  the  rate  of  change  of  the  dependent 
variables  Integrated  over  the  region  In  terms  of  their  values  on  the 
boundary  of  the  region  and  other  quantities  which  are  specified  for  the 
problem. 

A  typical  example  Is  the  one-phase  Stefan  problem.  Choosing  units  so 
that  the  latent  heat  Is  1  and  the  melting  temperature  Is  0,  one  may 


formulate  the  problem  classtcally  as  the  solution  of 


u^  ■  au  ,  X  c  n{t)  ,  t  >  0  ,  (l.la) 

JHTtFsSTtr  '  .  .  C-'b) 

u(x,t)  -0  ,  X  ^  R(t)  ,  t  >  0  ,  (1.1c) 

u(x,t)  UoCx)  ,  X  £  ,  (l.ld) 

where  u  1s  the  energy  per  unit  volume  and 

uo(x)  ^  1  ,  X  €  n(0)  ,  (1 .le) 

uo(x)  -  0  ,  X  ^  n{0)  .  (1  .If) 


To  complete  the  statement  of  the  problem,  one  must  Indicate  how  the  region 
n(t)  changes  with  time.  The  rate  of  change  of  n  Is  determined  by  requiring 
conservation  of  energy  and  noting  that,  with  the  units  chosen,  the  flux  of 

dll 

energy  Into  an  from  the  interior  of  n  Is  -  and  this  must  equal  the  flux 
of  material  volume  which  has  changed  phase  times  the  latent  heat.  Thus,  If 


V  Is  the  outward  normal  velocity  of  the  boundary  an(t), 

v»-|^  .  (l.lg) 

Since  conservation  of  energy  was  Invoked  to  complete  the  formulation 
of  problem  (1.1),  there  Is  a  certain  logic  to  writing  the  problem  directly 


as  a  conservation  law: 


u^  »  Af(u)  , 
f(u)  ■  max(u  -  1 ,0) 


(1.2a) 

(1.2b) 


u(t)  uo  .  (1.2c) 

Using  Gauss's  theorem,  we  see  from  (1.1c)  that.  If  n  Is  bounded,  |  udx  Is 
conserved.  It  Is  easy  to  check  that  the  left-hand  side  of  (1.2a)  has  a 
singular  part  and  the  right-hand  side  has  a  singular  part  - 
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where  is  a  Dtrac  measure  on  an.  Equating  the  two  singular  parts,  we 
recover  (l.lg). 

An  advantage  of  the  formulation  (1.2)  over  (1.1)  Is  that  In  (1.2)  the 
topology  of  the  set  {x|u(x)  ^1)  does  not  appear,  whereas  (l.lg)  Is  mean¬ 
ingless  If  n  Is  not  defined  on  3n(t)  or  Is  not  bounded.  A  more  subtle 
advantage  Is  that,  by  formulating  (1.2)  as  the  set  of  equations  (1.1),  we 
are  already  making  an  assumption  about  the  nature  of  the  solution,  and 
such  an  assumption  may  be  unwarranted.  In  (1.1),  we  assume  that  u  has  a 
jump  from  1  to  0  at  the  free  boundary,  that  Is,  that  the  front  between 
{x|u(x)  >^1}  and  {x|u(x)  =  0}  Is  sharp.  However,  In  the  Inhomogeneous 
problem 

u^  ■  Af(u)  +  a(x)  (1 .3a) 

where  o(x)  ^  0  Is  bounded  and 

supp  o  n  supp  Uq  0  ,  (1.3b) 

there  will  not  be  a  sharp  front.  Instead,  for  a  period  of  time  there  will 
be  a  diffuse  boundary  between  {x|u(x,t)  ^  1}  and  {xlu(x,t)  »  0),  and  any 
attempt  to  state  the  problem  In  the  form  (1.1)  will  be  doomed  to  failure. 

A  different  example  Is  afforded  by  the  hyperbolic  conservation  law 

u^  +  V  •  (f(u))  ■  0  ,  (1 .4a) 

u(x,t)  Uo(x)  .  (1.4b) 

Generally,  solutions  of  this  equation  will  not  be  continuous.  For  discon¬ 
tinuous  solutions  the  singular  parts  of  the  terms  on  the  left-hand  side  of 
(1.4)  will  have  to  cancel  one  another.  Define 

±  u*(x,t)  5  11m  supy^  ±  u(y,t)  (1.5) 

and 

8n(t)  2  {xlu'*’(x,t)  -  u*(x,t)  >  0>  .  (1.6) 
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Then  the  balancing  of  singularities  gives  for  the  velocity  v  of  3ft(t), 
v{u'*’(x,t)  -  u"(x,t))  •  f(u'*‘(x,t))  -  f(u"{x,t)).  (1.7) 

A  given  differential  equation  corresponds  to  a  multitude  of 
conservation  laws,  and  the  solution  of  the  equation  which  Is  considered 
acceptable  Is  critically  dependent  on  the  quantity  assumed  to  be  conserved. 
For  example,  the  equation 


may  be  written  as 


“t  ■*'  ““x  " 


(1.8a) 

(1.8b) 


or 

(u2)^  +  (|u3)  -0  .  (1.8c) 

In  (1.8b)  the  conserved  quantity  Is  |  udx.  In  (1.8c)  It  Is  |  u^x.  The 
general  solutions  of  the  two  conservation  laws  are  quite  different. 


2.  Numerical  Solution  of  Some  Simple  Conservation  Laws 


One  naturally  looks  for  algorithms  which  rigorously  preserve  the 
required  Invariants  of  the  conservation  laws,  and  In  addition  have  the 
feature  that  they  generate  sufficiently  smooth  solutions  of  the  conserve- 
tion  laws  with  satisfactory  accuracy.  For  the  hyperbolic  conservation  law 
(1.4),  such  an  algorithm  has  been  studied  (Ref.  5).  To  explain  how  the 
algorithm  works,  we  Introduce  the  sets 

£„(<*)  =  <x|u(x)  r  «}  ,  «>0  ,  (2.1a) 

Ey(a)  a  {x|u(x)  <  «)  ,  o  <  0  .  (2.1b) 

In  addition,  for  1  ^  1  £  N,  we  define  the  operators  S^(h)  by 

(S^(h)u)(x)  -  u(x  -  he^)  ,  (2.2) 
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» 


where  e^  Is  a  unit  vector  In  the  direction  of  the  x^-ax1s.  The  algorithm 
generates  a  collection  of  functions  u"(x),  n  ■  0,  1,  2,  •••,  dependent  on 
a  parameter  t.  called  the  time  step.  The  purpose  Is  to  approximate  u(x,t) 
by  u"(x)  for  t  ■  t/n  sufficiently  small.  u”(x)  Is  given  by 

u''(x)  =  G%o(x)  (2.3a) 

where 


Gu 


f  N 

n  S>(f'(a)  •  e.T)x(E  (a))da 
U  1=1  ^  1  U 

C  0  N 

n  S^(f'(o)  •  e^T)x(Ey(o))da 

J-OB  i  =  l 


(2.3b) 


Here  x(E)  Is  the  characteristic  function  of  the  set  E.  For  a  problem  1n 
one  space  dimension,  under  mild  regularity  assumptions  on  uo  and  some 
assumptions  on  f  (e.g.,  f  convex),  the  algorithm  (2.3)  has  been  shown  to 
generate  an  approximate  solution  with  an  error  0(t)  (Ref.  5). 

A  numerical  approach  to  a  more  general  class  of  conservation  laws 
which  Includes  (1.2)  and  the  one-dimensional  version  of  (1.4)  has  been 
presented  (Ref.  3).  One  considers  the  problem 

u^  +  Lf(u)  =  0  ,  (2.4a) 

u(x,t)  uo(x)  ,  (2.4b) 

where 

0  <  f(u)  -  f(v)  iu-v1fu>v  (2.4c) 

and  L  Is  a  linear  operator  such  that  the  semi-group 

S(h)  •  e“^^  Is  contractive  In  and  l"  (2.4d) 

The  numerical  method  generates  functions  u''(x),  supposed  to 
approximate  u(x,nT),  according  to  the  prescription 


u"{x)  ■  F%o(x)  .  (2.5a) 

Fu  5  u  -  f(u)  +  S(T)f(u)  .  (2.5b) 

The  functions  u”(x)  have  been  shown  to  converge  to  u(x,t)  as  ^  ■  t/n  -*■  0 
(Ref.  3).  For  the  special  case  of  problem  (1.2)  In  one  space  dimension, 
the  error  In  the  approximate  solution,  under  m11d  assumptions  regard¬ 
ing  the  regularity  of  uq,  has  been  shown  to  be  0^  In  (Ref.  2). 

Comparison  with  the  procedure  (2.3)  suggests  that  there  might  be  some 
theoretical  advantage  1n  solving  (2.4)  by  considering  the  algorithm 

u"(x)  -  G%o(x)  (2.6a) 

where 

Gu(x)  =  [  S(f'(o)T)x(E  (o))da  -  [  S(f ' (o)t)x(E„(o) )da 

io  “  “  (2.6b) 

At  this  point,  the  convergence  of  u”(x)  In  (2.6a)  to  u(x,t)  as  t  =  t/n  -►  0 

Is  only  conjectured.  It  may  be  of  Interest  to  note  that,  upon  use  of  the 

Interpolation  formula 

S(f‘(o)T)  1  +  f'(a)(S(T)  -  1)  ,  (2.7) 

(2.6)  goes  over  to  (2.5). 

3.  Function  Spaces  Approximate  to  Free  Boundary  Problems 

A  conservation  law  may  strongly  suggest  a  natural  function  space  In 
which  to  search  for  a  solution  of  the  associated  free  boundary  problem. 

This  Is  especially  true  when  the  conserved  quantity  can  be  thought  of  as 
a  norm  for  a  large  class  of  solutions.  Tie  space  has  already  figured 
In  error  estimates  for  approximate  solutions  of  (1.2)  and  (1.4). 

In  particular,  for  the  conservation  law  (2.4a),  It  follows  from  (2.4c) 
and  (2.4d)  that  the  operator  F  Is  contractive  In  L*.  Further,  the 
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algorithm  (2.5)  Is  stable  In  both  and  L*”.  In  cases  of  Importance,  F 
Is  stable  In  still  other  function  spaces.  For  example.  If  S(h)  not  only 
satisfies  (2.4d),  but  has  the  property  that 

(S(h)u)(x)  -  I  G(x' ;h)u(x  +  x' )dx'  ,  (3.1) 

one  may  find  a  whole  class  of  function  spaces  In  which  F  Is  stable.  One 
of  the  most  Interesting  Is  the  space  with  the  norm 

Hull  -  [  **  |3Ey(o)|da  .  (3.2) 

Using  this  space  and  the  algorithm  (2.5),  one  may  deduce  Immediately 
bounds  on  the  measures  of  surfaces  of  discontinuity  of  solutions  of  Stefan 
problems  and  solutions  of  model  hyperbolic  free  boundary  problems  In  terms 
of  the  Initial  data.  More  generally,  one  may  Introduce  function  spaces 
whose  norms  are  constructed  In  terms  of  capacities  of  sets  Ey(“).  Such 
spaces  and  others,  and  their  use  In  studying  the  regularity  of  solutions  of 
the  conservation  law  (2.4),  are  treated  In  more  detail  elsewhere  (Ref.  4). 

4.  The  Formal  Extension  to  Systems  of  Conservation  Laws 

One  might  consider  u  and  f(u)  to  be  n-vectors .  In  which  case  (2.4) 
becomes  a  system  of  n  conservation  laws.  The  numerical  scheme  (2.5)  still 
makes  sense  for  this  case,  and  It  Is  at  least  consistent  with  the  differ¬ 
ential  equation.  If  L  has  the  property  (2.4d),  and  If  In  addition 
Vu  e  r"  we  have 

|u  -  f(u)lj  +  lf(u)|j  ■  |u|j  ,  (4.1) 

then  we  get  stability  of  the  operator  F  In  L'; 

IIuIIli  ■  J  Iu|jdx  .  (4.2) 
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However,  the  contractiveness  of  F  In  L*,  which  followed  so  easily  for  the 
case  of  a  single  equation,  does  not  necessarily  follow  for  a  system  of 
equations.  In  addition  It  appears  that,  of  all  the  spaces  alluded  to 
above,  with  norms  constructed  In  terms  of  capacities  of  the  Ey(“).  the 
only  space  In  which  we  can  expect  to  get  ||F||  £  1,  except  for  trivial 
cases.  Is  Li. 

A  second  example  Is  afforded  by  Burgers'  equation  In  N  dimensions,  by 
which  we  mean  the  laws  of  mass  and  momentum  conservation  of  a  perfectly 
compressible  Inviscid  fluid  (no  pressure).  Burgers'  equation  Is 


+  V  ♦  (pu)  =  0  , 

(4.3a) 

(pu)^  +  V  •  (puu)  “  0  , 

(4.3b) 

p(x,t)  Po(x)  , 

(4.3c) 

(pu)(x,t)  (pu)o(x) 

(4.3d) 

A  formal  algorithm  seeks  the  solution  of  (4.3)  by  means  of  a  "Boltzmann 
equation"  approach  (Refs.  6,  8).  The  method  Is  similar  to  the  method  of 
(2.3).  One  writes 


(4.4a) 


where  G  Is  defined  by 


C) 


I  S^(v  •  e^T)|(p«(v  -  u))dv 
J  v^^n^  S^(v  .  e^T)j(p6(v  -  u))dv| 


(4.4b) 


and  1$  given  In  (2.2).  One  can  show  that  solutions  of  (4.3)  for  which 
u  Is  sufficiently  smooth,  uniformly  In  time,  are  given  exactly  by  (4.4). 
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However,  from  the  point  of  view  of  free  boundary  problems.  It  Is  the 
nonsraooth  solutions  which  are  of  primary  Interest. 

A  further  formal  extension  of  the  methods  described  above  would  be  to 
constrained  systems  of  conservation  laws.  By  these,  we  mean  simply  con¬ 
servation  laws  for  which  the  solution  Is  required  to  lie  In  some  subset  of 
a  function  space,  with  the  corresponding  d1fferentl!»l  equations  holding 
for  solutions  In  the  Interior  of  that  subset.  For  .xample,  the  equations 
of  Inviscid  hydrodynamics  may  be  described  succinctly  as  the  solution  of 
Burgers'  equation  (4.3)  subject  to  the  one-sided  density  constraint  p  £l 
(Ref.  6).  The  corresponding  free  boundary  problem  1s  the  water  wave 
problem. 

In  a  stepwise  solution  of  a  constrained  conservation  law,  It  Is  most 
natural  to  let  each  step  consist  of  two  procedures:  First,  the  solution 
for  a  step  of  the  corresponding  unconstrained  conservation  law,  as 
described  heretofore,  and  then  the  satisfaction  of  the  constraint.  In  the 
satisfaction  of  the  constraint,  any  numerical  procedure  should  carefully 
conserve  the  required  Invariants  of  the  system.  For  example,  in  hydro¬ 
dynamics  the  requirements  imposed  by  the  one-sided  constraint  on  the 
density  are  met  by  solving  another  conservation  law.  In  this  case  the  one- 
phase  Stefan  problem  (Refs.  6,  8). 

Another  example  of  a  constrained  conservation  law  Is  the  "diffusion- 
consumption"  problem  (Ref.  1).  Here  one  solves  the  conservation  law 

u^  "  au  -  1  ,  (4.5a) 

u(x,t)  uo(x)  ,  (4.5b) 
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subject  to  the  constraint 


u(x,t)  ^  0  .  (4.5c) 

In  some  cases  ft  may  be  fruitful  to  consider  the  specification  of 

N 

conditions  on  a  function  at  the  boundary  of  a  finite  region  ^  r  as  the 
Imposition  of  a  constraint  on  a  conservation  law  which  Is  formulated 
throughout  R^,  and  to  solve  the  resulting  constrained  conservation  law 
stepwise.  In  accordance  with  the  procedure  outlined  above. 

5.  Problems  Associated  with  Coupled  Equations  In 
Several  Space  Dimensions 

The  free  boundary  problems  of  mechanics  which  have  elicited  the 
greatest  Interest  have  Involved  coupled  equations  In  several  space  dimen¬ 
sions.  The  simplest  of  these  Is  the  water  wave  problem  described  above. 

A  slightly  more  difficult  problem  's  the  transonic  flow  problem,  which  can 
be  considered  to  be  given  by  the  conservation  laws 

+  V  •  (pu)  =  0  ,  {5.1a) 

(pu)^  +  V  •  (puu)  =  -  vp{p)  ,  (5.1b) 

where  P(p)  Is  specified.  The  free  boundary  here  is  called  the  sonic 
surface.  A  more  elaborate  system  of  conservation  laws.  In  which  another 
dependent  variable  (entropy)  and  an  extra  conservation  law  (for  energy)  are 
added,  describes  the  flow  of  an  Inviscid  compressible  fluid,  and  the 
corresponding  free  boundaries  are  shock  waves. 

For  each  of  these  sets  of  conservation  laws,  which  may  be  thought  of 
as  progressively  enriched  versions  of  Burgers’  equation,  similar  questions 
as  to  the  well-posedness  of  the  Initial  value  problem  arise.  We  noted  In 
the  last  section  that,  although  we  can  write  down  formal  algorithms  for 
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"solving"  these  conservation  laws,  the  variety  of  spaces  In  which  we  get 
stability  of  the  algorithms,  and  thus  strong  regularity  results,  diminishes 
drastically  on  passing  from  a  single  conservation  law  to  a  system  of  con¬ 
servation  laws.  That  the  proof  of  strong  regularity  results  for  such 
systems  may  not  be  just  difficult,  but  in  fact  Impossible,  is  strongly 
suggested  by  the  irregularities  actually  observed  in  the  physical  world, 
and  associated  with  the  names  of  Helmholtz  instability,  Taylor  instability, 
fingering  of  salt  water-fresh  water  interfaces  and  oil -water  interfaces, 
and  turbulence.  All  these  phenomena  are  associated  with  systems  of  equa¬ 
tions  in  more  than  one  space  dimension.  Since  Burgers'  equation  is  at  the 
core  of  all  these  sets  of  conservation  laws,  it  may  not  be  unreasonable  to 
expect  limitations  on  the  regularity  of  solutions  of  Burgers'  equation  to 
have  some  relevance  to  the  solutions  of  the  other  systems.  In  the  case  of 
Burgers'  equation,  it  would  appear  that  the  initial  value  problem  is  not 
well  posed  In  the  classical  sense  of  Hadamard,  and  that  it  may  be  neces¬ 
sary  to  go  to  a  concept  of  well-posing  in  a  stochastic  sense.  In  this 
case  It  may  be  that  a  given  initial  state  will  evolve  stochastically  after 
a  finite  time  (Ref.  7), 

Other  systems  of  conservation  laws  in  higher  dimensions  which  are  not 
built  around  Burgers'  equation  also  appear  to  have  highly  irregular  solu¬ 
tions.  As  a  prototypical  example,  we  might  consider  the  system 

u^  =  Af(u)  .  (5.2) 

For  example,  equations  like  (5.2)  have  been  used  to  describe  phase  changes 
In  materials  where  the  melting  temperature  depends  on  the  concentration  of 
a  solute.  In  some  of  the  models  proposed,  u  has  two  components. 
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representing  concentration  of  solute  and  energy  density.  It  would 
appear  that  the  free  surfaces  separating  phases  are  highly  unstable,  and 
that  this  Instability  Is  reflected  In  the  formation  of  highly  Irregular 
boundaries  containing  dendrites  (Ref.  9). 
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